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Abstract 

To describe the slow dynamics of a system out of equilibrium, but close to a dy- 
namical arrest, we generalize the ideas of previous work to the case where time- 
translational invariance is broken. We introduce a model of the dynamics that is 
reasonably general, and show how all of the unknown parameters of this model may 
be related to the observables or to averages of the noise. One result is a generalisa- 
tion of the Fluctuation Dissipation Theorem of type two (FDT2), and the method is 
thereby freed from this constraint. Significantly, a systematic means of implement- 
ing the theory to higher order is outlined. We propose the simplest possible closure 
of these generalized equations, following the same type of approximations that have 
been long known for the equilibrium case of Mode Coupling Theory (MCT). Natu- 
rally, equilibrium MCT equations are found as a limit of this generalized formalism. 

Key words: Fluctuation-Dissipation Theorem, glasses. Mode Coupling Theory, 
aging 
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Introduction 



Glasses are typically amorplious, dynamically slowed or arrested states of mat- 
ter, whose dynamical relaxational processes are dramatically slower than those 
of liquid states, despite their great structural similarities. It has long been 
known that, for many systems, slow and careful approach to what is clearly a 
non-equilibrium glass transition leads to entirely reproducible behaviour more 
reminiscent of phase transitions at equilibrium. Such phenomena are often 
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called, albeit loosely, 'equilibrium glass transitions'. The most complete cur- 
rent microscopic theory available for these systems, Mode Coupling Theory 
(MCT) [1] has implicit within it [2] the Fluctuation Dissipation Theorem of 
type two (FDT2) [3], but despite this, the agreement between theory and light 
scattering measurements [4] are remarkably good for simple colloidal systems. 
The explicit manifestation of the 'equilibrium' nature of the system is that all 
obscrvablcs arc functions only of the time differences, and are independent of 
when the experimental measurements are commenced. However, cases such as 
this are by no means natural, and typically experiments are found to depend 
on the time of waiting ty^ after a change in the external parameters of the 
system. Numerous experimental descriptions exist. However, in a pioneering 
paper [5] on model systems it was shown that equilibrium MCT ideas could 
be extended to deal with systems that 'age' and therefore violate FDT2 [6-8]. 
The ideas contained in that paper have been helpful in developing the con- 
cept of aging, but it is not until recently that attempts have been made to 
extend the existing microscopic theory of equilibrium glassy states (MCT) in 
the same spirit. Two such approaches appear to be developing. One of these 
is described in this paper. The other [9] uses an entirely different strategy by 
which projection operators are extended to the non-equilibrium state. 

Before beginning our discussion, we should sound some notes of caution. 
Whilst the framework we build is quite general, the approximations we use are 
the analogues of those used in equilibrium MCT [1], and in non-microscopic 
theories of aging [5], [10]-[13]. We have argued elsewhere [2,14] that these 
amount to a type of dynamical approximation, and that is known to be ac- 
curate close to the arrest transition only for colloidal systems. Elsewhere, as 
in molecular glasses, they are approximately correct up to some characteristic 
distance from the glass transition, often called the MCT temperature [15]. In 
terms of approximations, what we here present cannot be expected to improve 
on this fundamental limitation. On the other hand we may hope that what is 
described below will be the non-equilibrium (aging) theory of those systems 
for which MCT has proven of value. 



1 Fundamentals 

It has been shown in previous works [2,14] that it is possible to write density 
equations of motion starting from Newtons equations. This gives an explicit 
formula for Pk(^) where the density variables Pk{t) are defined as the Fourier 
transform of the number density p{r,t), i.e. 




(1) 
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In principle, to describe the underlying Newtonian variables of position and ve- 
locity wc would need many fields, the most fundamental of which are generally 
considered to be the hydrodynamic variables of longitudinal and transverse 
currents, along with a local temperature or entropy variable. Until now, mode 
coupling type theories have been based on longitudinal current, and therefore 
only the density degree of freedom, and we shall illustrate our points in the fol- 
lowing within this particular approximation. Thus, the first time derivative of 
density is related to the longitudinal current, and the derivative of this current 
is treated in Mode Coupling Theory (MCT), leading to an equation involving 
the second order derivative of density in time [1]. Extensions involving more 
fields are more complex, and they will be the subject of future work. 

The basic approach in any case is to develop a general model that can de- 
scribe slow out-of-equilibrium dynamics by extending FDT2 [3], and then to 
find some simple approximations to close the model, as a starting point of 
applicability. 

Thus, we write the density equations as. 



where we separate out deterministic and stochastic motions, representing with 
T>]i{t, t) and TZi^i, r) respectively the deterministic and the random force acting 
on the density waves. Observables have to be calculated by averaging over the 
noise distribution which has not yet been stipulated. 

We have explicitly indicated the dependence on two distinct times in the his- 
tory of the system, t and r. In general, t can represent the final observation 
time, while r can be considered the so-called waiting time in an aging ex- 
periment. The choice of these two times implies that, once the variables of 
the system are known at time r, equations (2) determine the dynamics of the 
density variables through the unknown functions 'Dk(^,'?') and TZk{t,T). 

We start by discussing some general properties that the random force must 
possess. Indeed, it is fundamental to require the stochastic process to satisfy 
the causality condition. Essentially all that we will show, up to the final ap- 
proximations, derive only from the form of equations (2), and the causality 
requirements. These explicitly are. 



where t corresponds to all times later than r. The brackets indicate averages 
over the (unknown) non-equilibrium distribution of the system. If the decom- 



Pk(i) = 2^k(^,T) + 7^k(^,T) 



(2) 



(p_k(r)7ek(t,r))=0 

(p-k(T)7^k(^,T)) = o 



(3) 
(4) 
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position into deterministic and stochastic motions, made in (2), was exact, 
these averages would be the same as the averages over the noise distribution. 



Now, using (2), the two causahty conditions (3) and (4) can be rewritten as, 

(p_k(r)pk(t)) = (p-k(r)I)k(t,r)) (5) 
(p_k(r)pk(t)) = (p-k(r)Pk(t,r)) (6) 



These also imply two further conditions, that are obtained by taking their first 
derivative with respect to r, and then using (3) and (4) as well as (2). Thus, 
we have, 

(p_k(r)^Pk(i,r))=0 (7) 



d 

(p_k(r)— I?k(t, r)) = (p-k(r)7^k(^, r)) 

= {V_^{t, T)n^{t, r)) + (7^_k(r, r)7ek(t, r)) (8) 



Conditions (7,8) constitute two of the fundamental constraints on which to 
build a theory of non-equilibrium slowed dynamics. 

At this point we are still at liberty to use any trial deterministic force. We 
choose to write the most general linear form, i.e. 

V^(t, t) = -Qk WPk W - Mt)pi.(t) 

t t 

- I dt'^i,{t,t')pi,{t') - J dt'xi.{t,t')f^{t') (9) 

T T 

where the quantities f2k(t), t'k(^),7k(^5 and Xk(^,^') are four unknown pa- 
rameters of the theory. We may note that so far, if one accepts the underlying 
Mori-type hypothesis that dynamics may be represented by a deterministic 
(slow) part, and noise (fast) part, even to describe non-equilibrium states, the 
linear approximation could be essentially exact for most conditions, providing 
the memory kernels are chosen correctly. 

Substituting (9) into equations (2), we obtain the set of generalised Langevin 
equations for the density variables as, 

Pk(t) + i^k(t)Pk(t) + ^^k(t)Pk(t) 
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t t 

T T 



(10) 



With the expression (9), we assume that the deterministic force is composed 
of hnear terms in pk(^) and pk(^), exphcitly separated into instantaneous 
and non-local contributions, l^k(^) and ^'k(t) are time dependent because the 
system is evolving towards equilibrium and, therefore, so are the collective 
variables. f2k(t) is the frequency of the elementary excitations that would 
be 'phonons' at equilibrium, while vi^it) takes into account the dumping of 
the modes. The integral terms are memory contributions to these kinds of 
effects, configurational evolution and dissipation respectively. They are not 
time-translational invariant as the observables of the system are not. 

It is not possible to apply a simple FDT2 as soon as this invariance is broken. 
Nevertheless, the choice made in (9) enforces the width of the noise distribution 
to be closely related to the memory kernels as was for the equilibrium case. 
This can be viewed as a generalization of FDT2 for a system out of equilibrium. 
Indeed, combining (7,8) and (9) we obtain. 



(7e_k(r, r)7^k(^, t)) = n{ J,(r) - --^ 7k(t, r) (11) 




and 



Xk(t,T) = -]p^{t,r)^\ogSk{r) (12) 



with the definitions. 



^fc(^,T) = ^(p_k(r)pk(t)) (13) 

Jk{t,T)^^{p.^{T)p^{t)) (14) 
Su{t) ^ Sk{T,T); Jk{r) ^ Jk(r,r) (15) 



Here, S'fc(r) and Jfc(T) are equal time correlators. In the equilibrium case, 
clearly, they simply reduce to, 

St\t)=S,; Jt'\r) = ^ (16) 

pm 



Also,5^)(t) = 0. 
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Thus, equation (11) is the generahsation of FDT2 for systems out of equi- 
hbrium, and indeed it possesses the correct hmit to the equihbrium FDT2 

[3], 

(7^Lt(r,r)7^(-)(t,r)) = ^lt'\t - r) (17) 



while Xk*^(i, t) = 0, and this quantity does not appear in the equihbrium 
theory. 

It is still necessary to determine the generalised frequency, f2k(t), and the 
instantaneous part of the friction, ^'k(i), to characterise completely the model. 
Thus, we write equations (10) for the time r, 

Pk(r) + z/k(r)pk(r) + l^k(r)pk(T) = n^{T, r) (18) 



and we require the following initial conditions for our stochastic process, 

(p_k(T)7^k(T,T))=0 (19) 

(p_k(r)7^k(r,r)) = 0. (20) 

These are consistent with the causality relations of equations (3,4), but one 
can imagine other models where they are not applied. Their choice leads to 
the following relations. 



hk{r) - Mr) + \Mr)Sk{T) + f^k(T)5fe(T) = (21) 

iMr) + Mr)Jk{r) + In^Skir) = (22) 
Solutions to these equations are, 

' ^Jk{t)S,{T) - SI{t) ^ ' 

= , r r NO / X — 7;T^-^ (24) 



^Jk{r)Sk{T) - Sl( 



Note that the choice (9) implies four unknown parameters, and equations 
(7,8,23,24) fix these in terms of observables, and of the noise distribution, 
both of which might be considered input into the theory or determined self- 
consistently. It is interesting to note that a very similar structure has been 
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developed by Latz [9], using the method of non-equihbrium projector opera- 
tors. 

We shall discuss the first attempts to determine these quantities self-consistently 
within certain approximations. 

Following the same kind of steps made in [2,14] we define, 

Mt) + Mt)Mt) + ^i^{t)Pk{t) = Mt) (25) 



where ^k;(^) are considered to be the couplings between the modes. If we seek 
to minimize these, then the 'best' choice of modes can be determined by. 



^ "^m') = (26) 



^ (l^k(t)r) = (27) 



These equations give the same solutions as the ones given in (23) and (24), 
providing one means of interpreting the choice (19) and (20). 

Now we define the normalised correlators. 



, . Sk{t,T) (P-k(r)pk(t)) . . 

and, using the Langevin equations (10) we write. 



r^2 Q 

-Q^^kit, r) + Ui,{t)—^k{t, t) + Ok(t)$fc(t, r) 

+ 1 dt'x^it, t')^k{r, t') + j dt'^^it, t')g^Mr, t') = (29) 



^2 Q 

* * d 

+ / dt'x^{t,t')^k{r,t') + j dt'^^{t,t')—^k{r,t') = (30) 

r T 

(29) and (30) are the equations of motion for the system. This set of equations 
is equivalent to the equations proposed in [9]. 
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Now, these equations require knowledge only of the dispersion of the noise 
{TZ{t,T)TZ{T,T)) (see equations (11,12)), as well as Sk{T) and Jkir). 



Before turning to make useful approximations, we discuss briefly a new method 
that leads to closure of these equations. It can lead to a sequence of corrections, 
and, in the equilibrium case leads to the known MCT equations within the 
same type of approximations. 

The idea consists in applying a variational principle for that part of the exact 
random force that is hard to calculate. We shall write a sort of constitutive 
relation for the noise that connects it to density variables, but where the noise 
is constrained to have certain reasonable properties, essentially ensuring that 
both sides of the stochastic equations are consistent. We therefore write, 

T^Vit.r) = nu{t)pu{t) + i/k(t)pk(t) + d^^\t)pu{t) + ^ d£U^)Pk-k'(^)Pk'(t) 

k'^k 

+ 4'!k',k" (^)pk-k" (t)Pk"-k' (t)Pk' {t) + ... 

k',k"^k' 
t t 

+ j dt'^^{t,t')p^{t') + J dt'xk{t,t')p^{t') (31) 



T 



where the set of parameters have to be determined by making this trial 
form the closest possible to the true one, that contains p]s_{t). In principle, one 
needs an infinite number of the terms in the sum (31), and the level of the 
approximation of a theory will be correspondent with the number of terms 
considered. 

To be acceptable this choice of noise must be 'faithful', in the sense that the 
Langevin process of equations (10) should be preserved. This issue can be 
partially addressed by requiring that the noise fluctuations on a single time 
slice are consistent with the Langevin process (10). Thus, we impose, 

^ "7^k(t,r)-7^r^«(^,r)H = (32) 



dd(^\t) 

where TZk{t, r) is defined in (10). This leads to an infinite set of coupled equal- 
time equations. In the next section we shall give an example of these. 

The other conditions to be satisfied by the expansion (31) of the noise are the 
constraints (19,20). These conditions can now be explicitly written as. 



f^k(r)5,(r) + lMr)Sk{r) + di\T)S^{T) + ^ ^ 41(T)(p-k(T)pk-k'(T)Pk'(r)) 

k'j^k 



8 



+^ Yl 4%,k"(^)(P-k(r)pk-k"(^)Pk"-k'(r)pk'(r)) + ... = (33) 
\^^.{T)S,{T) + v^{t)J,{t) + \di\T)S^{T) + 1 ^ <U^)(P-k(r)pk-k'(r)Pk'(r)) 

+ ^ E 41,k"(T)(P-k(T)pk-k"(t)pk"-k'(T)pk'(T)) + ...=0 (34) 

These would be exact equations of motion determining the non-equihbrium 
structure factor and currents, providing one knows the couphngs S"'^ and the 
expansion (31) is exact. Wc now discuss some simple closure ideas. These 
should not be regarded as complete, but merely indicative of the sorts of 
approximations that will need to be considered in future. 



2 Simplest Closure 



It is now necessary to employ some approximations regarding the multiple av- 
erages involving densities, or equivalently the noise distribution, as well as the 
kinetic terms. These last terms would be well approximated by their equilib- 
rium limits if we are exploring the system after the velocities have equilibrated, 
and this is the line of thinking we shall pursue as a first approximation. A dif- 
ferent approach has to be taken to account for the multiple averages and also, 
for the equal time structure factor that is itself evolving in time. 

To make relation to the known MCT theory at equilibrium we shall consider 
only the first non-linear term in the sum which represents the noise (31), and 
thus, the only parameters left to choose are d^\t) and d^\^,{t). We then apply 
the minimization condition (32), 



d 



ddi\t) 
d 



|pk(t) - di\t)p^{t) - E 41(^)Pk-k'(^)Pk'(t)r ) = (35) 

k'jtk I 



p^{t) - di\t)p^{t) - E <Ui)pk-k'(^)Pk'(^)r ) = (36) 



^C?k^p(^) \ k'j^k 



This explicitly gives the conditions, 

E<L'W(p-k(t)Pk-k'(^)Pk'(^)) 

k'jtk 
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(37) 



E4,ic'W(P-k+p(i)p-p(t)Pk-k'(t)Pk'(t)) 

+ 4'^(t)(p_k+p(t)p-p(t)pk(t)) - (p-k+p(t)p-p(t)pk(t)) = (38) 

to solve for finding the best possible coefficients d^\t) and d^^i^,{t). Equation 
(37) is trivial, as it corresponds to the same equation as (33). 

If we had chosen to work in the Newtons equation representation we would at 
this point still have the bare potential present in the problem [2,14]. In essence, 
by choosing the best form of the noise in terms of density in an independent 
manner, as above, we renormahse the instantaneous forces experienced by the 
density waves at a single time slice, away from what they would have been 
if we had the bare potential to generate the noise correlations [2,14]. This 
may also be viewed as performing a partial sum over some of the noise to 
give exact, or nearly so, equal-time averages, and leaving the averages to be 
calculated in the remaining correlators to be carried out with the remaining 
noise. The most important approximations are then due to approximations of 
the correlations between different time slices. In a manner, we may view this 
as having renormalised the vertices of the problem prior to proceeding with 
any approximation of the unequal time correlators. 

At this stage we could in principle solve for the coefficients, but consistent 
with previous MCT approximations, we may further simplify the problem. As 
in equilibrium MCT [1] , we choose a Gaussian approximation for the four point 
averages and a superposition approximation for the triplets, (p_k(i)Pp(^)Pq(^)) ~ 

p+q, which consists in neglecting the triple direct correla- 
tion function C3. 

Recall that these averages are not over all the noise ensemble but only that 
part remaining after the pre-averaging mentioned above. In this sense, using 
the leading non-trivial approximation for averages (Gaussian for even averages, 
superposition for odd averages) may not be as severe an approximation as one 
would think at first sight. 

Indeed, in the equilibrium case solutions of equation (38), in these approxi- 
mations, are given by, 

where cj^'^^ is the equilibrium direct correlation function. We have also assumed 
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j(2){e9) q{eq) q{eq) qieq) _ 
l^k'+k "k,k' '-'-k '-'k-k"^k' — 

Thus, inserting this formula in (31) we can calculate (11) neglecting the inte- 
gral contribution, but still accepting (29) as the correct equations of motions, 
recovering the well established equihbrium MCT equations [1,14]. 

In principle, out of equihbrium, it is stiU possible to solve for the coefficients 
(ik k,(t), that minimize (36), but to have some explicit expressions for them we 
have to make further approximations than those in the equilibrium case. We 
shall, as a first approximation, consider ourselves to be in the regime where 
the kinetic contributions to the observables are the same as in the equilibrium 
limit, and we thereby assume that the velocities relax towards their equilibrium 
values much faster than the positions. 

In the following we consider the case of any time instant t', in the range 
between r and t. Thus, we approximate. 



(p-k(t')Pk-p(t')Pp(t')) « - ^^'^^ ^^ ^p(0 - ^^'^II'-pKO- (40) 



The last approximation to be made is that {p-k{t')pk{t')) — {l/2)Sk{t') — 
Jk{t') « -k'^/pm, with Sk{t') ^ 0. 

We can now write the result for the best coefficients in this approximations. 



with nckit') = 1 — 1/Sk{t') the generalised direct correlation function, with n 
the number density N/V. Evidently, this is a sort of adiabatic approximation, 
in which, ultimately, we expect changes in the noise distribution to arise from 
changes in the slowly evolving structure. 

We now use (10) and (31) to obtain. 



^k(i', r) ^ Yl (k • k')cfe'(OPk-k'(i')Pk'(i') + Mt)pUt) 

t t 
+ J dt'j^{t,t')p^{t') + J dt'x^.{t,t')p^{t') (42) 

r T 

where t' can be any instant between r and t. The term containing the memory 
kernel is neglected in equihbrium MCT [1,14]. We will, therefore, also neglect 
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the two integral terms, which in any case would be of higher order in our 
arguments. 



Now, if we approximate the four-point density correlation for different times, 
by Gaussian decomposition, as before, and we neglect the contributions arising 
from the term z/k(t)pk(^) in (42) because of higher order, we finally find the 
generalised expression for the memory kernel 7k(^, t). 



and the memory kernel Xk(^,T) is then obtained combining this expression 
with (12). We can now apply successive levels of approximation to equation 
(33, 34) to find an equation of motion that determines the non-equilibrium 
structure factor, S'k(t). In this way, we have a simple theory of the non- 
equilibrium dynamical structure factor. 



Conclusions 

In this paper we have derived a generalised Langevin equation that should 
reasonably describe an out-of-equilibrium dynamically slowed system. This 
leads to new constraints on the noise distributions rather than traditional 
FDT2, but these constraints are still practicable to apply. We then expand 
the noise in a 'Landau type' expansion of the density variables, leading to a 
non-linear generalised Langevin process. However we conceive that some pre- 
averaging of the variables has taken place so that exact equal time averages 
are recovered to a high level of quality by assuming remaining noise is only 
gaussian distributed. To preserve fidelity of the Langevin process we then 
insist that the noise be a faithful representation, and that the structure of the 
Langevin process is still preserved. These requirements are implemented by a 
variational principle to determine the renormalised coefficients of the noise, 
and by implementing the orthogonality of the noise to the 'slow' variables, 
Pk{t) and Pk(i)- This leads in principle to the determination of all of the 
unknown coefficients of the noise. 

Further approximations are possible but not obligatory, but they make the 
contact with existing ideas easier. The outcome is that we end up with a the- 
ory that must be considered as an 'adiabatic' extension of the existing MCT 
in the sense that kinetic terms are neglected, assumed to be equilibrated. The 
MCT equations are therefore modified, but the form of the memory kernel 
is preserved, only it now includes the dependence on the waiting time. In 



7k(^, r) 




-F(k- k')k- (k- k')cjk/(T)c|k-k'|(^)]'S'|k-k'|(^,T)5'jk'(i,T). 



(43) 
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eqTiilibriTim MCT one typically inputs any good approximation to the static 
structure factor for which the YBG (Yvon-Born- Green Equation) equation of 
constraint should hold, or hold approximately [2]. Here we deduce the corre- 
sponding conditions, and these generate an equation of motion for the equal 
time structure factor and current correlator. These equations would be exact 
providing the expansion chosen for the random force, as in (31), is a complete 
one. As non-equilibrium structure factors to use as inputs are much less known 
than equilibrium ones, these equations, or some similar ones, will be crucial 
to close the theory. 

The strategy we have pursued here seems different to that developed in [9], 
and it is as yet too early to see how these different approaches will relate to 
each other. Probably the answer to that question is that in the end, as with 
the equihbrium MCT case [14], a fair degree of correspondence will emerge, 
though it will be important to see if either approach leads naturally to higher 
levels of approximation. 

Prom what we know about colloidal systems near their kinetic arrest, this the- 
ory should be reasonably successful in describing the phenomena qualitatively, 
perhaps, as with equilibrium MCT, quantitatively. This opens the possibility 
to begin systematic study of the aging of colloids and soft materials [16,17]. 
However, to make useful progress, in the regime where FDT2 is violated for 
molecular glasses, it will require better approximations than those shown here, 
and this is a matter to which we, and others, most certainly will direct our 
attentions. 
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